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O ; ABSTRACT 

^ ■ Aims. As part of the Herschel guaranteed time key program 'HOBYS', we present the PACS and SPIRE photometric survey of the star forming 

^ [ region Vela-C, one of the nearest sites of low-to-high-mass star formation in the Galactic plane. Our main objectives are to take a census of the 

„ cold sources and to derive their mass distribution down to a few solar masses. 

■I Methods. Vela-C has been observed with PACS and SPIRE in parallel mode at five wavelengths between 70 yum and 500 yum over an area of about 
3 square degrees. A photometric catalogue has been extracted from the detections in each of the five bands, using a threshold of 5 cr over the local 

' background. Out of this catalogue we have selected a robust sub-sample of 268 sources, of which ~75% are cloud clumps (diameter between 

. 0.05 pc and 0.13 pc) and 25% are cores (diameter between 0.025 pc and 0.05 pc). Their Spectral Energy Distributions (SEDs) have been fitted 

Mh ■ with a modified black body function. We classify 48 sources as protostellar, based on their detection at lOfim or at shorther wavelengths, and 218 
as starless, because of non-detections at 70 fim. For two further sources, we do not provide a secure classification, but suggest they are Class 

f-i ' protostars. 

^ ^ Results. From SED fitting we have derived key physical parameters (i.e. mass, temperature, bolometric luminosity). Protostellar sources are 

I . in general warmer ({r>=12.8K) and more compact ({diameter>=0.040pc) than starless sources ({r>=10.3K, {diameter>=0.067 pc). Both these 

Q ■ evidences can be ascribed to the presence of an internal source(s) of moderate heating, which also causes a temperature gradient and hence a more 

" peaked intensity distribution. Moreover, the reduced dimensions of protostellar sources may indicate that they will not fragment further A virial 

, analysis of the starless sources gives an upper limit of 90% for the sources gravitationally bound and therefore prestellar in nature. A luminosity 

■ vs. mass diagram of the two populations shows that protostellar sources are in the early accretion phase, while prestellar sources populate a region 

' of the diagram where mass accretion has not started yet. We fit a power law A'(logM) oc M"' '*"-^ to the linear portion of the mass distribution of 
prestellar sources. This is in between that typical of CO clumps and those of cores in nearby star-forming regions. We interpret this as a result of 

" ' , the inhomogeneity of our sample, which is composed of comparable fractions of clumps and cores. 
^ , 

, Key words. ISM: individual objects: Vela-C - ISM: clouds - stars: formation - submillimeter: ISM - (stars): circumstellar matter 

^: 

I • 1. Introduction (MS X, IRAS) and around thirty isolated protostars and seven 

' , ^ , . ^ , ^ . , , , embedded young clusters associated with C'^O clumps have 

Most of the stellar content m our Galaxy forms in ^^^^ ^^^^^ ^Liseau et al. 1992; Lorenzetti et al. 1993; Massi 

(M . 30 K) and dense (« > 10^ cm- ) cores withm the Giant Molecular 2003, Baba et al. 2006). A bright HII region, RCW 36, 

^ . Clouds (GMC) of the Galactic Disk (e.g.. Blitz 1991). The main ^^^^ associated by Massi et al. (2003) with an early-type 

. ; properties of the stellar population, such as the efficiency of the ^^^^ ^^^^^1 Q5_gQ^ ^ ^^^^^^ .^^^ 34^ 

. mass conversion into stars and the shape of the initial mass func- originally associated to Vela-C by several authors (Herbst 1975, 

;XJ ■ tion (IMF) may very well be closely related to the physical prop- ^ ^ ^ ^991), but recent observations favour a much 

H . erties and mass distabution of the progenitOT structures m the ^.^^^^^^ ^^^2.5 kpc, Bik et al. 2010). Vela-C has been 

. . parental cloud (e.g. Mc Kee & Ostriker 2007). 250yum, 350^m and 500 ;um with the Balloon-borne 

Being one of the nearest GMCs in the Galactic disk the Vela Large Aperture Submillimeter Telescope (BLAST, Pascale et al. 

Molecular Ridge (VMR, Murphy & May 1991, / . 260°-275% 2OO8) that provided the first census of tJie compact dust emission 

b^+3), represents an ideal observational tax-get. It is composed ^ evolutionary stages and lifetimes (Netterfield et al. 

of four molecular clouds at a distance between 700 pc (clouds 2009) 
A, C, D) and 2000 pc (cloud B, Liseau et al. 1992), and harbours 

protostars with masses up to 10 M^, both isolated and clustered ^ target of the Herschel guaranteed time key program 

(Massi et al. 2000, 2003). Vela-C is the most massive compo- "OBYS ( Herschel imaging survey of OB Young Stellar ob- 

nent and hosts the youngest stellai- population (Yamaguchi et J^'^^s', Motte et al. 2010), Vela-C has been observed with the 

al. 1999); this latter has been investigated by combining neai-- (Poglitsch et al. 2010) and SPIRE (Griffin et al. 2010) 

infrared data (J, H, K-hand images) with far-infrared surveys '^^"^^i'^^ between 70A<m and 500 A/m. The extended emission in 

the form of filaments and ridges has been presented by Hill et 

* Herschel is an ESA space observatory with science instruments pro- al. (201 1) who identify five diff^erent sub-regions with different 

vided by European-led Principal Investigator consortia and with impor- column densities. Here we will focus on the determination of 

tant participation from NASA. the physical parameters of the compact sources, their evolution- 
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ary stage and their mass distribution. The latter has been inves- 
tigated in several star forming regions through different tracers 
(e.g. Testi & Sargent 1998, Motte et al. 1998, 2001, Kramer et 
al. 1998, Reid & Wilson 2006, Enoch et al. 2008, Rathborne et 
al. 2009, Konyves et al. 2010, Ikeda & Kitamura 2011). In this 
paper, we intend to enlarge the statistics providing the mass dis- 
tribution over a wide range, going from subsolar values up to 
tens of solar masses. 

This paper is organized as follows. In Section|2]we describe 
the observations, the data reduction procedures and the pho- 
tometric results. In Section[3]we describe the Spectral Energy 
Distribution (SED) analysis along with the fitted physical pa- 
rameters. In Section|4]we discuss the evolutionary stages of the 
sources along with their mass distribution. Our results are then 
summarized in Section|5] 



2. Observations, data reduction and results 

Vela-C was observed on May 18, 2010 with the Herschel par- 
allel mode, i.e., using simultaneously PACS at 70/160 jum and 
SPIRE at 250/350/500 yum. A common area of ~ 3 square de- 
grees was covered by both the instruments around the position 
a(y2000.0)=08h59™55\ V2000.0)=-43°53'00". The field was ob- 
served in two orthogonal directions at the scan speed of 20"/s. 
The data reduction strategy is described in detail in Traficante 
et al. (2011): here we summarize only the fundamental steps. 
From archival data to the Level 1 stage, we used scripts prepared 
in the Herschel Interactive Processing Environment (HIPE, Ott 
2010), partially customized compared to the standard pipeline. 
The obtained time ordered data (TODs) of each bolometric de- 
tector were then processed further by means of dedicated IDL 
routines and finally maps were created using the FORTRAN 
code ROMAGAL. Images of Vela-C at all wavelengths have 
been presented by Hill et al. (2011), here for the readers' con- 
venience, we show in Figure [T| the composite 3 -color image at 
70/160/500 /zm. 

2. 1 . Source Detection and Photometry 

The detection and photometry of compact sources have been 
carried out using the Curvature Threshold Extractor package 
(CuTeX, Molinari et al. 2011) applying the same strategy as 
Molinari et al. (2010). Briefly, in each band, sources were iden- 
tified as peaks in the second-derivative images of the original 
HOBYS maps, then an elliptical Gaussian fit was performed to 
provide: 1) the total flux, integrated down to the zero intensity 
level, 2) the observed FWHM (defined as the geometrical mean 
of the major and minor ellipse axes), and 3) the peak intensity. 
The latter, when divided by the local rms noise, allows us to ob- 
tain an a posteriori estimate of the S /N ratio. We filtered out all 
sources with S /N < 5. Following Elia et al. (2010), entries at 
different wavelengths in the PACS/SPIRE catalogue have been 
attributed to the same source based on simple positional crite- 
ria. In practice, we associate two sources detected in two differ- 
ent bands if their mutual angular distance does not exceed the 
radius of the Herschel half-power beam-width (HPBWQ) at the 
longer wavelength. Around 15% of the entries in the catalogue 
present multiple associations at decreasing wavelength; as a gen- 
eral rule, the closest counterpart was kept without attempting to 
divide the flux at the longer wavelength among the sources. We 



' The values of the Herschel HPBW are 5.0" at lOfim, 11.4" at 
160 Afm, 17.8" at 250/im, 25.0" at 350 yum, and 35.7" at 500 /zm. 



assigned to each source the celestial coordinates of the counter- 
part at the shortest wavelength, where the spatial resolution is 
higher Finally, to exclude artifacts, we discarded sources with 
an axes ratio > 2 and a position angle randomly changing with 
wavelength by more than 20°. 

Table 1. Statistics of the 5 cr PACS/SPIRE catalogue. 



Band # Sensitivity limit 90% Completeness limit' 



(Jy) (Jy) 

Total entries" 1686 - 

70yum 658 0.04 0.21 

160 yum 871 0.09 0.67 

250yLim 966 0.11 1.05 (17) 

350 yum 697 0.35 1.32(22) 

500 yum 416 046 1.95 



" For comparison, in parentheses the BLAST completeness limit are 
reported, estimated from Netterfield et al. 2009 (their Figure 3). 
Number of sources detected in at least one band. 



2.2. Catalogue statistics 

The complete Herschel catalogue of Vela-C will be published 
in the next future as part of the HOBYS data products. Here 
we limit to give in Table [T] a statistical summary of the detected 
sources. In the second column we list the number of entries in 
each of the five bands: this number increases with increasing 
wavelength up to /I = 250 yum, than decreases at longer wave- 
lengths. Such behaviour can be explained in the light of the 
sensitivity limits given in the third column: since the number 
of entries in the first three bands increases despite the loss of 
sensitivity, such increasing is a real effect reflecting the intrinsic 
properties of the cloud. Conversely, in the last two SPIRE bands, 
the poorer sensitivity limits, together with the loss of angular 
resolution, are likely responsible for the reduced number of the 
detected sources. 

The flux completeness limits (fourth column) are estimated 
at the 90% of confidence level, estimated by recovering artifi- 
cial Gaussian sources randomly spread over the map and whose 
fluxes and diameters cover the same ranges found for the ob- 
served sources. In Figure|2l we show the percentage of the recov- 
ered artificial sources as a function of the flux for each Herschel 
band. For comparison, the BLAST completeness limits at 250 //m 
and 350/im (estimated in Figure 3 of Netterfield et al. 2009) are 
more than a factor of 15 larger than the SPIRE limits at the same 
wavelengths. 

2.3. Source diameters 

In each of the five bands, the physical FWUMdeci^l) has been 
derived by deconvolving the observed FWHM with the HPBW 
at the same wavelength: 

FWHMdeciX) = -yjFWHMl - HPEWj. (1) 

Noticeably, V^llMdeciX) increases with A, with mean values 
of 7.5", 13.4", 22.7", 28.6", 41.0" in the five Herschel bands, 
respectively. Since we are most interested in tracing the cold dust 
emission (namely that at temperature ^ 25 K), which is likely 
poorly related to the 70yum flux, we define as angular source 
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Fig.l. Composite 3-color image of Vela-C : PACS 10 /im (blue), PACS 160 A^m (green), and SPIRE 500 jum (red). The blueish 
regions near the center and towards the bottom right of the map are the HII regions RCW 36 and RCW 34, respectively. Red and 
blue diamonds represent the locations of prestellar and protostellar sources, respectively, belonging to the sample of 268 objects 
used for the SED analyisis. Bigger diamonds indicate sources with M > 20 M©, while yellow diamonds indicate a few sources inside 
RCW 34 (see text). The top left segment corresponds to a length of 5 pc in the sky for a distance to Vela-C of 700 pc. 
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Fig. 2. Percentage of recovered artificial sources as a function of 
their fluxes in the five Herschel bands. 



diameter (6) the FWUMdec measured at 160 yum, where we detect 
a considerable number of sources with a good angular resolution 
(see also Motte et al. 2010). For a small fraction of sources that 
remained undetected at 160//m, 9 is the FWHM^/fc measured at 
250^m. 

We consider as spatially resolved the detections with 6 larger 
than ~6Q% of the angular resolution at 160jum (correspond- 
ing to a diameter D ~ 0.025 pc at the Vela-C distance), as also 
adopted in the Aquila Rift and Polaris by Andre et al. (2010). 
Noticeably, around ~ 25% of our sources have a diameter be- 
tween 0.025 pc and 0.05 pc ('cores'), while most sources are 
larger ('clumps', 0.5 < D < 0.13 pc); in the following analysis, 
we will not distinguish any longer between the two categories, 
and will refer to both of them with the general term 'source' . 



3. Analysis 

3.1. Flux scaling 

The physical parameters of the detected sources have been de- 
rived by fitting their SEDs with a modified black body curve. 
Two examples of the observed SEDs are shown in Figure[3] 
These illustrate a trend seen in many of the SEDs, namely that 
the observed fluxes (depicted with crosses) flatten or even rise 
longwards A > 250 jum. This effect is a direct consequence of 
the increase of the FWHMded^) with wavelength (Sect. l2.2] i. 
which in turn implies increasing areas over which the emis- 
sion is integrated. In practice, such flattening indicates that a 
single-temperature, modified black body fit is not completely 
adequate to model the observed photometric points. However, 
since the small amount of available data prevents us from using a 
multiple-temperature, modified black body model, we restricted 
the fit only to the innermost portions of the sources, whose spa- 
tial scales are defined by the source angular diameter 0. To es- 
timate the emission coming from these restricted solid angles, 
we followed the flux scaling method adopted by Motte et al. 
(2010) and described in detail by Nguyen Lifo'ng et al. (2011). 
This method is based on the idea that for quasi-spherical, self- 
gravitating sources, the radial density law is oc r^^ (with r ~ 0.1-1 
pc), and thus M(< r) oc r (Shu 1977). This relationship implies 
that fluxes scale linearly with the source radius, provided that the 
emission is optically thin and the temperature gradient is weak 
(T(r) oc r"*, with q ~ Q, see Motte & Andre 2001 for more de- 
tails). All these conditions have been empirically verified in our 
case by taking the photometry at different apertures for some 
isolated sources, and finding a linear increase of fluxes with in- 
creasing radii. Hence, fluxes at 250/vm, 350//m, and 500 jum, 
were reduced by the ratio FWYiMdecWIO- The average values 
of these factors are 1.9 at 250;um, 2.0 at 350//m, and 3.2 at 
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Fig. 3. Examples of modified black body fits to the SEDs of 
two sources (undetected and detected at IQfim, respectively). 
Crosses and dots mark the data points before and after the cor- 
rection applied at longer wavelengths to minimize the effects of 
the poorer spatial resolution (see text). Values of the main fitted 
parameters are reported. 

500 //m, respectively. In Figure|3] dots represent fluxes after the 
above correction. 



where Fy is the observed flux at the frequency v, Q is the 
source area (in sr), 'QyiJd) is the black body function at the 
dust temperature Td, and Ty is the optical depth, parametrized 
as Ty - (v/vo'f, where vq - c/Aq is the frequency at which 
T=l, and /3 - 2, as predicted by simple dust emission models 
(Hildebrand 1983). A linear least-squares fit to the scaled fluxes 
(with their photometric uncertainties) was performed by com- 
paring the observations with a database of models where Tj and 
/!() are allowed to vary in the following ranges: 8 K < < 40 K; 
lOyum < Aq < 500 jum. We applied the PACS color corrections, 
amounting to a maximum of 4% of the flux, to the model fluxes 
before comparing them with the data (Miiller et. al 2011). As a 
check, we also verified that the fitted values of Ao never exceed 
160 //m (see Sect. l4.1l and Figure|5]i; this insures that the emis- 
sion is optically thin in the SPIRE wavelength range, namely 
where flux scaling is performed (see the discussion in Sect. 13. II ) 
Finally, masses were derived from the dust emission using 
the relation 

M = (t/2n/^300)T300 (3) 

where ^300= 0.1 cm^ g"' is the opacity per unit gas mass 
computed at A = 300 //m assuming a gas-to-dust ratio of 100 
(Beckwith et al. 1990), and d-100 pc. (The details of the deriva- 
tion of Eq.[3]from Eq.|2]will be given in Pezzuto et al., in prepa- 
ration). 



3.2. Selection criteria 

Once the SPIRE fluxes were rescaled, we further filtered the cat- 
alogue to select a robust sub-sample for SED analysis. The final 
list was composed of sources fulfilling the following criteria, i.e., 
sources: 

(1) with detections in at least three adjacent bands between 
160jum and 500 /zm, to contain the peak of the emission of 
the cold dust; 

(2) without a dip in the SED between three adjacent wave- 
lengths; 

(3) not peaking at 500 yum. Note that conditions (2) and (3) ex- 
clude most of the artifacts or detections affected by incorrect 
positional associations or bad photometry; 

(4) spatially resolved at 160/im; 

(5) not presenting multiple associations at A> 160yL(m, to allevi- 
ate confusion problems; 

(6) not belonging to the RCW34 region, to exclude sources 
that could be more distant than 700 pc (see Sect.[T]l. These 
sources, identified by eye in the Herschel maps, are depicted 
in yellow in Figure[T] 

The first three criteria are the most stringent ones, since they 
lead to a selection of only 388 sources out of 1686 present in the 
catalogue. Applying criteria (4) to (6), we get a final sample of 
268 sources to be fitted. Their locations are depicted in Figure 
[U with different colors indicating their evolutionary stage (see 
below). 

3.3. SED fitting 

Since cold dust emission is likely unrelated to the photometry 
measured at 70 fim, we fitted the SEDs longward /I > 160 jum and 
checked a posteriori whether or not the modified black body 
curve fits also the 70 fim point. The assumed relation is: 

Fy = (l-e-^OBv(W (2) 



3.4. Separating starless from protostellar sources 

Since Vela-C was not observed by Spitzer at 24 yum, which would 
have revealed embedded sources with adequate sensitivity, we 
used the Herschel data themselves, and where available, other 
mid- and far-IR surveys to distinguish starless and protostellar 
sources. This distinction was primarily done on the basis of the 
presence of a 70 yum counterpart (see also Bontemps et al. 2010). 
Indeed, a tight correlation between 70 //m fluxes and internal lu- 
minosities of protostars has been demonstrated (Dunham et al. 
2008), and, furthermore, external heating by the interstellar ra- 
diation field cannot produce a 7Qfj.m emission from a source 
detectable with Herschel. Hence, we consider as starless all 
sources without a 70;um counterpart (e.g. Figure[3] left panel). 
Conversely, protostellar are: /) sources whose 70//m flux lies 
above the best modified black body fit by more than 3 cr (e.g. 
Figure[3] right panel), or, //) sources detected in the mid-IR atA< 
70 /im by other surveys (IRAS/MSX/Akari). Out of 268 objects, 
218 are starless and 48 are protostellar Two further sources, not 
detected at A < 70 /im by other surveys, were fitted with a single 
modified black body curve from 70 fim to 500 yum. The temper- 
ature of both sources is ~ 25 K, i.e. warmer than starless sources 
(Sect. l3.4T l. On this basis we surmise that these sources are em- 
bedded Class protostars (Andre et al. 1993). 

3.5. Physical parameters 

In this section, we discuss the statistics of the key physical 
parameters measured directly from the observations or derived 
from SEDs fits. These quantities are listed in Tableland plotted 
in Figure|4] 

Starless sources exhibit a very narrow temperature distribu- 
tion (top left panel), strongly peaked at ~ 10 K. The high 
cutoff is at Td ~ 15K, while the low one coincides with the 
minimum temperature of 8 K considered in the fit. On aver- 
age, protostellar sources are warmer than starless sources; their 
temperature distribution spreads almost uniformly between 9 K 
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and 15 K, though a few sources are found with up to 25 K. 
The overall temperature distribution (not shown in Figure|4|l re- 
sembles that found from BLAST data of the whole Vela region 
(Netterfield et al. 2009). This latter distribution, however, does 
not show the sharp ~ 10 K peak, but rather has a shallower shape 
peaking at ~ 12 K. This marginal difference may be due to both 
the better angular resolution of Herschel, which allows deeper 
penetration into colder regions, and to the inclusion of warmer 
regions (e.g., Vela-D, Olmi et al. 2009) in the Netterfield et al. 
sample. It is also interesting to compare the compact source dis- 
tribution with the average temperature of the surrounding cloud. 
The latter was derived by Hill et al. (201 1) by fitting a modified 
black body to the Herschel SED of each pixel. In the parts of the 
cloud where no compact sources are detected, the temperature is 
fairly constant at ~ 14 K (black line in Figure|4|. Practically 
all the starless sources are colder than the ambient medium, a 
result expected for dense structures heated only by external ra- 
diation. Interestingly, the same result is found for most of the 
protostellar sources. This finding reasonably indicates that the 
external envelope we are tracing is moderately heated by the in- 
ternal protostar(s), which is therefore likely to be in the very 
early phases of the evolution. 

The mass distribution (top right panel of Figure|4]i shows a 
nearly monothonic decrease from lower to higher masses. There 
are no significant differences between the average masses of star- 
less and protostellar sources (~ 5 M©, Table|2|; this suggests that 
the envelope remains substantially unperturbed by the internal 
gravitational collapse and, again, that the protostellar sources 
are relatively unevolved, as little mass accretion is taking place. 
Noticeably, 35 sources (of which 7 protostellar and 28 starless) 
have 10 < M/Mq < 20 and 8 (of which 1 protostellar and 7 star- 
less) have 20 < M/Mq <55, whose locations are depicted with 
big crosses in Figure[T] Apart from a few exceptions, these mas- 
sive objects are rather compact (with diameters between 0.04 pc 
and 0.07 pc), hence they are potentially able to form individu- 
ally massive stars. For comparison. Hill et al. (2011) found in 
Vela-C nine objects with masses between 13 Mq and 70 Mq. We 
consider the agreement between the two determinations satis- 
factory, given the different source extraction methods (Hill et al. 
2011 takes into account in the photometry the bias created by 
bright rimmed effects and source multiplicity) and the different 
prescriptions under which the SED fitting has been performed 
(for example. Hill et al. did not apply flux scaling). 

The source diameter distributions are presented in the bot- 
tom left panel of Figure|4] For starless sources, it spreads from 
the resolution limit of 0.025 pc up to 0.13 pc, peaking around 
0.05 pc. Such dimensions, as anticipated in Sect.|2l span from 
those of progenitors of single stars, up to those typical of large 
clumps, that may give rise to the formation of protostellar clus- 
ters. Conversely, protostellar sources are significantly more com- 
pact; this suggests that no further fragmentation will occur and 
that the majority of the protostellar sources likely host single 
stellar objects or stellar systems with a small number of mem- 
bers. The small sizes of protostellar sources can be alterna- 
tively (and additionally) explained by considering that the in- 
ternal source(s) of heating likey causes a temperature gradient 
and hence a more peaked intensity distribution that is fitted with 
a compact spatial extent. 

Finally, in the bottom right panel of Figure|4l we plot the 
distributions of luminosity in the Herschel bands (Lfir)- This 
quantity has been computed by interpolating the corrected fluxes 
between all available wavelengths between 70 jum and 500 yum. 
For starless sources, Lfir is a good approximation of bolometric 
luminosity; typically, it ranges between 0.05 L© and 0.5 L©, with 




D (pc) Log Lp,u(L3) 

Fig. 4. Histograms of the main physical parameters of the 
sources in Vela-C (top left: temperature; top right: mass; bot- 
tom left: diameter; bottom right: Herschel luminosity). Starless 
and protostellar sources are separately highlighted in red and 
blue, respectively. The black line in the top left panel represents 
the temperature of the surrounding cloud, taken from Hill et al. 
(201 1), while that in the bottom left panel highlights the spatial 
resolution limit. 

a small percentage of objects at larger luminosities not exceeding 
5 Lq. For protostellar sources, LpiR has to be considered as a 
lower limit to the bolometric luminosity because of the lack of 
near- and mid-infrared data points in the SEDs. 

4. Discussion 

4.1. Virial analysis of starless sources 

Starless sources are defined as prestellar if they are gravitation- 
ally bound (e.g., Andre et al. 2000, Di Francesco et al. 2007), 
hence they can potentially form one or more stars. In principle, 
spectroscopic observations would be required in order to derive 
the virial masses and to securely establish the sources' dynami- 
cal states. Since such observations are not available for Vela-C, 
we have assumed thermal pressure support and neglected the in- 
ternal turbulence. In this simplified view, the virial mass can be 
surrogated by the critical Bonnor-Ebert mass: 

Mbe ^ l.AR^^a^lG (4) 

Here a is the sound speed at the source temperaturfl G is the 
gravitational constant, and Kbe the Bonnor-Ebert radius (in pc). 
We take Rbe - Djl. Starless sources with M/Mbe > 0.5 are se- 
lected as being gravitationally bound (Pound & Blitz 1993) and 
candidate prestellar. By this definition, a large number of prestel- 
lar sources is found (206 objects), i.e., ~ 94% of the starless 

" a= -^JkeTj/^, where ks is the Boltzmann constant, is the fitted 
temperature, and jx = 2.33 m(H) is the mean molecular weight, with 
m(H) the atomic hydrogen mass. 
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Table 2. . Statistics of the main physical parameters of the sources in Vela-C. We give median, average, minimum and maximum 
values. 



Parameter 




Starless (218) 




Protostellar (48) 




median 


average 


min-max 


median 


average 


min-max 


M(Mo) 


3.3 


5.5 


0.13-55.8 


2.7 


4.8 


0.15-29.1 




10.0 


10.3 


8.0-15.2 


11.4 


12.8 


9.0-24.2 


Z)(pc) 


0.064 


0.067 


0.025-0.13 


0.040 


0.040 


0.025-0.07 


i'FIR (Lq) 


0.17 


0.22 


0.04-4.8 


0.6 


8.0 


0.08-138 
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Fig. 5. Mass vs. diameter diagram for the sources in Vela-C. 
Red filled circles and red open circles are prestellar and unbound 
sources, respectively. Blue filled circles are protostellar sources. 
Dotted lines are the loci of sources with M = 0.5 Mbe, for T = 8 
K and T - 20 K, while the dashed vertical line represents the 
spatial resolution limit. We also show the curve representing the 
relation between mass and diameter in the case A()-160 fim. This 
relation shifts towards higher masses with increasing Aq: hence, 
since all our sources are located below the reported curve, they 
have optically thin emission at A> 160 yum. 

ones, remarkably higher than the percentage (69%) found in 
the Aquila Rift cloud (Konyves et al. 2010, Andre et al. 2010). 
Moreover, this percentage mantains high (193 objects, i.e., 88%) 
even if we demand that bound sources have M at least equal to 
IMbe- 

To check further this result, we plot all sources in a 
mass vs. diameter diagram (Figure|5]l. First we note that, as a 
consequence of condition (4) in the selection criteria of Sect. 13. 21 
all points lie to the right of the spatial resolution limit of 0.025 
pc, shown with a vertical dashed line. Second, we plot the re- 
lation between mass and diameter as a function of the sec- 
ond variable in the Equation[3] namely the wavelength Ao at 
which the optical depth t-1. This relation shifts towards higher 
masses with increasing Aq, therefore all the sources located be- 
low (above) the curve for a given value of Aq have optically thin 
(thick) emission at A > /Iq. As example, we show such relation 
for the value Aq =160/im, which is the shortest wavelength that 
we include in the SEDs fits. Since all our sources are located 
below this curve, their emission is optically thin at wavelengths 
longer than 160/zm. 



Consistently with our virial classification, most of the 
prestellar sources in Figure|5] lie above or between the lines of 
0.5 Mbe, computed at 8K and 20 K. We note, however, that a 
relevant fraction of these objects have a mass up to a factor of ten 
larger than the Bonnor-Ebert mass at 20 K. Reasonably, such ob- 
jects are dynamically unstable, if turbulence and magnetic field 
supports against gravity are neglected. 

All the sources classified as unbound are located below the 
Bonnor-Ebert curve at T = 8 K. Noticeably, most of them clus- 
ter around the value M ~ 0.2 Mq. This latter corresponds to the 
mass detection limit, estimated from the flux sensitivity limits of 
Table[T]and by assuming average values of T and Aq in Eqs.(2) 
and (3). Since such mass limit decreases with increasing tem- 
perature, we deduce that we are able to efficiently probe only the 
unbound sources particurlarly warm. Indeed, for the sub-sample 
of this category (12 objects) we get (rrf)=13.8 K, namely higher 
than the value of (r^)=10.3 K, referring to the whole sample of 
starless sources (see Table|2]l. Noticeably, our inability to detect 
cold, unbound sources, may explain the very large fraction of 
prestellar sources with respect to the unbound ones (see above). 
Moreover, this fraction could further decrease if the most mas- 
sive prestellar cores were gravitationally unstable, as already no- 
ticed above. 

Finally, to show more clearly the results of Table|2] we also 
plot in the same diagram the positions of protostellar sources. 
As already noted in Sect l3.5l they cluster in Figure|5] at diame- 
ters smaller than starless sources, and close to our spatial reso- 
lution limit. Since no significant differences in mass (as well as 
in temperature values) are recognizable with respect to starless 
sources, this result should reflect an increase in average density. 

4.2. Luminosity vs. Mass diagram 

Another interesting perspective to be examined is the relation 
between the bolometric luminosities (Lboi) and the envelope 
masses of the sources (M). To make such comparison, we take as 
Lb„i the Lfir values, even if these latter well represent bolomet- 
ric luminosities only for starless sources. Previous studies (e.g. 
Saraceno et al. 1996, Andre et al. 2000, Molinari et al. 2008, Elia 
et al. 2010, Bontemps et al. 2010, Henneman et al. 2010) have il- 
lustrated that the Li,oi vs. M plot is a meaningful tool for the char- 
acteiization of the evolutionary status of cores and clumps, since 
tracks representing model predictions can be directly compared 
with the locations of the observed sources. In such a diagram, a 
collapsing core is expected to follow initially an almost vertical 
path (i.e., mass constant with time, and luminosity strongly in- 
creasing with accretion), then subsequently protostellar outflow 
activity increasingly disperses the envelope while luminosity re- 
mains almost constant, resulting in a horizontal track. 

In Figure|6] the Lhoi, M pairs of the Vela-C protostellar and 
prestellar sources are plotted, along with the evolutionary tracks 
for low-mass objects adopted by Molinari et al. (2008), based on 
the simplified assumption of a central protostar accreting mass 
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Fig. 6. Lhoi - M diagram for prestellar and protostellar sources in 
Vela-C (red and blue circles, respectively). Grey solid lines rep- 
resent the evolutionary tracks for low-mass objects adopted by 
MoHnari et al. (2008), for initial values of 0.5, 1, 5, and 10 Mq, 
respectively. An arrow indicates the evolution direction, while 
green lines delimit the region of transition between Class and 
Class I sources (Andre et al. 2000). The red dotted line repre- 
sents the best-fitting power law (L^o/ °^ ^0-42±o.04-j ^^j. jj^g (jjs^-j. 
bution of the prestellar sources. 



from the envelope at a constant rate M = 10"^ yr '. More 
evidently than in previously shown diagrams, a significant dis- 
tinction between the two populations, although not a complete 
segregation, is seen. In the same plot, we delimit with green 
lines the transition region between Class and Class I objects 
(Andre et al. 2000), which roughly separates the vertical por- 
tions of the tracks, where Class objects are expected to be 
found, from the horizontal ones, where Class I objects are lo- 
cated. Noticeably, the large majority of the protostellar sources 
populate the region corresponding to Class objects. This is in 
contrast with our previous results based on SED fitting, accord- 
ing to which 48 sources are young protostars (although not more 
evolved) and just two sources are identified as candidate Class 
sources. Such a discrepancy can be however partly reconciled if 
considering that our underestimate of the bolometric luminosi- 
ties of protostellar sources reflects underestimates of the actual 
ages. Moreover, the Andre et al. (2000) limit separating Class 
and I could not be completely adequate in Vela, as already found 
in Taurus by Motte & Andre (2001). 

Also, the distribution of the prestellar sources is quite flat 
and homogeneously concentrated in a region of the diagram cor- 
responding to a very early stage where no mass accretion has 
started. A tentative power-law fit gives a dependence Lb„i oc 
^({).42±o.04) rpjjjg slope is shallower than that found by Brand et 
al. (2001) for their sample of protostellar sources, namely 0.85. 
A similar scaling relation has been found for the CO luminosity 
of clumps: from the theoretical point of view, Leo ^ M~} at the 
virial equilibrium (Scoville & Sanders 1987; see also Wolfire et 
al. 1993 and references therein); this slope has been found by 
Rengarajan (1984), while other authours find lower values (e.g. 
Solomon et al. 1987, Yonekura et al. 1997, Miyazaki & Tsuboi 
2000). A direct comparison with our result, however, cannot be 
performed due to the different classes of sources considered, or 
different quantities compared. 



the on-going star formation in Vela-C 




-0.5 0.0 0.5 1.0 1.5 2.0 

Log (M/Ma) 



Fig. 7. Source mass distribution in Vela-C. The error bars corre- 
spond to y/N statistical uncertainties. Continuous red line repre- 
sents the fit to the linear portion of the mass distribution (namely 
that consistent with a single straight line within the errorbars) 
for all the prestellar sources. A similar fit performed on the re- 
stricted sample of prestellar sources with diameters less than 
0.08 pc is represented with a grey dashed line. The derived slopes 
are reported as well. Mass distribution of CO clumps (Kramer et 
al. 1998), the single-star IMF (Kroupa 2001) and the multiple- 
system IMF (Chabrier 2005) are shown for comparison with an 
orange, a green, and a blue line, respectively. The dotted vertical 
line represents the mass completeness limit of prestellar sources. 



4.3. Source mass distribution 

FigureQ shows the mass distribution of prestellar sources in 
Vela-C. We estimate a mass completeness limit of ~ 4 M©, de- 
rived by assuming the average temperature of prestellar sources 
of Table|2] and the flux 90% completeness limit at 160 //m (see 
Table[T]). Note that this value is significantly better than the 
BLAST completeness limit (~ 14 M© for sources colder than 
14K). We fit a power law A^(logM) oc to the linear portion 
of the distribution, finding a slope y =-1.1+0.2. The uncertainty 
was determined by considering both the statistical error of the 
data and the variation of the slope with the histogram binning, 
which was varied from 0.1 to 0.3 in logiM/Mo). Noticeably, no 
variations are found in the slope taking into account only the 193 
objects for which M > IMbe, since their mass distribution dif- 
fers from that of the whole sample only in the mass range below 
the completeness limit, where the fit is not performed. 

The value of y we find in Vela-C is shallower than that found 
by Netterfield et al. (y=-1.9±0.2, for sources with Tj < 14 K). If 
we limit the fit to the BLAST mass completeness limit of 14 M©, 
however, we obtain a slope y = - 1 .9 ± 0.2, which reconciles our 
high-mass end distribution with that obtained with BLAST. More 
important, and thanks to the HerscheVs, massive increase in sen- 
sitivity, which significantly extends the mass range, a possible 
change of slope is recognizable at M ^ 10 M©. This is close to 
the value of ~ 9 M© where a steepening was probed, for exam- 
ple, in Orion A (Ikeda et al. 2007). Recent theoretical models by 
Padoan & Nordlund (2011) predict such change of slope, even 
if at lower masses (M ~ 3 - 5 M©). Unfortunately, no firm con- 
clusion can be drawn on the base of our data because of the poor 
statistics in the high-mass bins. 
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Compared with other literature values, our y value of -1.1 
is between the mass distribution slope of CO clumps (y — 0.7, 
Kramer et al. 1998) and that typical of prestellar cores, as for 
example, that measured in the Aquila Rift (y = -1.45±0.2) on a 
sample of sources with 2 ^ M/Mq $ 10 and diameter typically 
less than 0.08 pc (Konyves et al.2010). Hence, our slope re- 
flects the heterogenous nature of our sample, which is composed 
of objects with a variety of diameters typical of both clumps 
and cores. Indeed, if we include in the mass distribution fit only 
sources with diameter less than 0.08 pc (142 objects), we obtain 
a slope 7=-1.4±0.2 (grey dotted line in Figure|7]), which well 
reconciles with the slope of the Aquila Rift. Such a slope is also 
consistent with the IMF slope for 1.0:$ M/Mq (for single-star 
IMF y=-1.3±0.7, Kroupa 2001, while for multiple-systems IMF 
7=- 1.35+0. 3, Chabrier 2005), likely indicating that fragmenta- 
tion will proceed less efficiently in these small objects. 

Finally, a substantial agremeent within the error bars is 
provided by the comparison with surveys covering a mass range 
largerly overlapping that of Vela-C. For example, y ~ -1.3±0.2 
for 0.8 < M/Mq < 6 in Perseus, Serpens and Ophiucus (Enoch 
et al. 2008) and y = -1.3+0.1 for 3 ^ M/Mq ^ 60 in Orion A 
(Ikedaetal. 2007). 



5. Conclusions 

We have reported on the Herschel observations of the Vela-C 
star forming region over an area of ~ 3 squared degrees. From 
our analysis we obtained the following results: 



- A luminosity vs. mass diagram for the two populations 
of prestellar and protostellar sources has been constructed. 
Prestellar sources cluster in a well defined region of the dia- 
gram corresponding to a very early stage in which no mass 
accretion is expected to have started. A tentative power-law 
fit to the observed distribution gives a dependence Lhoi oc 
^(().42±o.04) Conversely, protostellar sources populate the di- 
agram region corresponding to the early accretion phase. 

- The mass distribution of the prestellar sources with M > 4 
Mq shows a slope of -1.1 ± 0.2. This is between that typi- 
cal of CO clumps and those of cores in closeby star-forming 
regions, maybe reflecting the heterogeneous nature of our 
sample, which is a mixture of cores and clumps. We signal a 
possible a change of slope in the mass distribution for M <: 
10 Mo, even if the very big errorbars in the higher mass bins 
prevent us to draw firm conclusions on its reliability. 
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- From a 5 cr level catalogue of cold compact sources in five 
Herschel bands between 70yum and 500 /im, we have se- 
lected a robust sub-sample of 268 sources. Their physical 
diameters indicate that our sample is mainly composed by 
cloud clumps (0.05 pc ^ D ^0.13 pc), together with a 
~ 25% of cores (0.025 pc <, D ^ 0.05 pc). 

- Based on the detection of a 70yum flux, we have identi- 
fied 218 starless and 48 protostellar sources. For two further 
sources we do not give a secure classification, but suggest 
them as candidate Class protostars. 

- Source physical parameters have been derived from modi- 
fied black body fits to the SEDs. Both starless and protostel- 
lar sources are on average colder than surrounding medium. 
This indicates that the radiation from interstellar field and/or 
from embedded protostars is unefficient in penetrating the 
cold dust in deep. 

- Protostellar sources are on average sligthly warmer and more 
compact than starless sources. Both these evidences can be 
ascribed to the presence of an internal source(s) of moderate 
heating, which also causes a temperature gradient and a more 
peaked intensity distribution. Moreover, the reduced dimen- 
sions of protostellar sources may indicate that they will not 
fragment further 

- No significant differences are found between the masses of 
the two groups. These range from sub-solar up to tens of so- 
lar masses. In particular, we find 8 objects with M > 20 M©, 
which are potential candidate progenitors of high-mass stars. 

- More than 90% of the starless sources result prestellar (i.e. 
bound) if a virial analysis is applied. This percentage, how- 
ever, should be considered as an upper limit, both because 
our sensitivity does not allow us to efficiently probe the cold- 
est unbound sources and because a number of the massive 
prestellar sources could be gravitationally unstable. 
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